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A new approach is proposed to the solution of the quantum mechanical inverse scattering problem 
at fixed energy. The method relates the fixed energy phase shifts to those arising in an auxiliary 
Sturm-Liouville problem via the interpolation theory of the Weyl-Titchmarsh m-function. Then a 
Marchenko equation is solved to obtain the potential. 
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INTRODUCTION 

Inverse scattering theory based on the spectral prob- 
lem of the s— wave Schrodinger equation has been worked 
out in the fifties [H-Q ^-nd developed further in the eight- 
ies with the aim to apply to problems in atomic, nuclear 
and subnuclear physics P, Q . This type of inverse scat- 
tering theory provides the fixed- Z potentials derived from 
input spectral data (eigenvalues and normalization con- 
stants) as well as scattering phase shifts 5i{k) known for 
all energy E — {hkY /(2m). 

Another type of inverse method has been developed in 
the sixties |6|-l8| to recover the fixed energy (or fixed-k) 
potentials from a set of phase shifts {(5i=o,i,...,imax} given 
at a particular wave number k. For practical analysis of 
the measured scattering angular distribution, the modi- 
fied Newton-Sabatier (mNS) method proved to be a pow- 
erful procedure Also the Cox- Thompson method, 
developed in the last decade to practical problems, offers 
results with nice physical properties such as non-singular 
behavior of potential at the origin or demanding only a 
finite set of input data [III [l3 to recover the potentials. 

In this letter a third type of inverse method will be pre- 
sented. It belongs to the fixed energy category requiring 
a set of phase shifts at one energy but uses spectral infor- 
mation belonging to an auxiliary problem resulted from a 
Liouville transformation. The new method bears resem- 
blance to a recent development by Horvath and Apagyi 



(HA) [13|, However, whereas HA obtain results from 



the Gelfand-Levitan (GL) equation by solving a moment 
problem, the present method makes use of the Marchenko 
integral equation in connection to the auxiliary spectral 
problem. 

Because in all physical problem the potential can be 
known beyond a finite radius r = a < oo it is enough for 
a practical inverse scattering method to seek the inverse 
potential q{r) within a finite range < r < a [beyond 
which q{r > a) = is assumed]. The HA method applies 
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a variable transformation r — aexp(— a;) and thereby 
casts the finite interval [0, a] in r to the half line [oo,0] 
in X. The inverse potential is obtained from the deriva- 
tive of the transformation kernel K{x, x) which is deter- 
mined from the GL equation within the interval [x,Q]. 
This means that the potential q{r') will be known in the 
interval [r, a] corresponding to the transformed interval 
[x, 0]. In order that the potential be known at the origin 
r = 0, the kernel K(x, x) should have been calculated 
accurately as x — >■ cx). This, however, is difficult in view 
of the restricted number of input data. 

If one tries, contrary to the HA method, to solve the in- 
verse problem by using the Marchenko integral equation 
then one gets the transformation kernel K{xtX) within 
the interval [oo, x\ which suggests that a more accurate 
potential value can be expected at the origin r — Q than 
in the case of the HA method. As we will see this task is 
feasible by straightforward numerical evaluation of cer- 
tain integrals. The basis of the theory remains the same 
as in the HA method, i.e. the recognition of the fact that 
the input set of fixed energy phase shifts {^i=o,i,...,;„ax} 
is related with the m-function of the auxiliary (Liouville- 
transformed) inverse spectral problem. While the HA 
method leads to the solution of a moment problem where 
the moments are calculated also from the input phase 
shifts, the present procedure requires the Jost functions 
to be determined for the auxiliary problem. Restricting 
ourselves to the simple case when there exists no bound 
state in the auxiliary spectral problem, we present some 
illuminating examples and point out the necessary future 
development related to the generic case. 

In section 2 the formalism will be outlined. Here an 
interesting (and hitherto seemingly overlooked) connec- 
tion between the GL and Marchenko equations is also 
discussed together with an adopted approximation of the 
m-function. Section 3 contains the illustrative examples, 
in a part where analytic solutions are known for the aux- 
iliary (transformed) problem and, in another part where 
there is no such a pre-known information available. Nev- 
ertheless, in all the cases we present also the results which 
are obtained without the use of the analytic solutions. 
These examples can help us to assess the accuracy of the 
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underlying numerical procedure and the approximations 
used. Finally, section 4 is devoted to the discussion and 
conclusion. 



THEORY 



Consider the radial Schrodinger equation (in units of 



h^l2m = 1) 



+ q{r) - 



ifi{r)=-l{l+l)ipi{r), re(0,a), 

(1) 



which, after a Liouville-transformation 
y(x,-a + l/2)2) = r-i/V/M, 



(2) 



takes the form of a standard Sturm-Lioville (or s— wave 
Schrodinger) eigenvalue equation 



?/(x, A) = Ay(x, A) xe(0,oo) (3) 



with the transformed potential 

Q{x)^r'^{q{r)-e), r 
The m-function of the Sturm-Liouville (SL) operator 
<P 



(4) 



Q{x) 



is known at the point values A = — (/ + 1/2)^, I = 
0, 1, 2, . . . , Zmax in terms of the phase shifts Si generated 
by the potential q{r) [with q{r > a) = 0] at a fixed wave 
number k, because the following relation [l^ holds: 



{ (' + 1/2) . 2)2) 



2/(0, -(? + 1/2)2) 

■//+i/2(fcQ) - tan 5iY{^^/^{ka) 
Ji+i/2{ka) - ta.ndiYi^i/2{ka) ' 



-ka 



(5) 



and implies the relation 



^ dJC{x, x) 
dx 



(9) 



between the potential and the transformation kernel 
which can be obtained upon insertion of ([7]) into Q. 

The input function J-{x) to the Marchenko equation 
© contains contributions of the bound states (discrete 
spectrum) and scattering states (continuous spectrum) 

^ 1 1 P'^ / \ 

(10) 

where B denotes the number of the bound states sup- 
ported by the potential Q{x), and Aj < 0, rrij > 
denote, respectively, the eigenvalues and normalization 
constants of the eigenfunctions y{x, Aj) belonging to the 
Schrodinger equation ([3]). The phase function A(k) cor- 
responds to the (s— wave) phase shift arising when a par- 
ticle with wave number k is scattered by the hypothetical 
potential Q{x). 

It is related to the modulus of the Jost function 
/+(k) — |/+(K)|e'^^''-' by the dispersion relation [l^ 



1 r 



\0g\ f+{K')\dK' 



K — K 



(11) 



The modulus, in turn, is connected to the ni- function by 



— lim 



1 



K > 0. (12) 



i.o+ ImTO(K2 + is) ' 

The latter equality can easily be proved by using the re- 
flection principle m[z) = m{z) and the Wronskian prop- 
erty f~{f^y — if^Yf^ ~ ^.ik of the two linearly inde- 
pendent Jost functions 

In summary, if the m-function is known near the posi- 
tive half axis of its argument one can determine the po- 
tential Q (or q). 



Connection to the Regge-Loeffel-Sabatier-Levitan 
framework 



The Marchenko equation reads as 



J'{x+y)+IC{x,y)+ K.{x,t)F{t+y)dt = Q {x < y) 

Jx 

(6) 

for the transformation kernel )C{x,y) from which the 
wave function can be obtained by the Povzncr-Levitan 
representation as 

1 1 f°° 

y{x, \) = -j= iiin{V\x) + ^(^' y) sin{V\y)dy. 

(7) 

This assumes a Dirichlet boundary condition 

y(0,A) = 0, 2/'(0,A) = l, (8) 



Let us transform back both the input and transforma- 
tion kernels into the r-space. Performing the transforma- 
tion X — — log ^ and using the notations 



/C (r, r' ) = IC ( — log — . — log — 
a a 



(13) 



J-(r, /) = ^ log ^ - log y = ^ log ^ j (14) 

we get from the Marchenko integral equation © the fol- 
lowing equation of Gclfand-Levitan type 

r ~ ~ dr" 

T{r, r')+/C(r, r')+ / /C(r, r")^(r", r') =0 (r > /). 

Jo r" 

(15) 
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The (wanted) fixed energy potential can also be calcu- 
lated as 



q{r) 



2 d/C(r, r) 
r dr 



(16) 



Note that equations (IT5|) and (fTBl) are exactly the same 
as those appearing in the Regge-Loeffel-Sabatier-Levitan 
framework of inverse scattering theory at fixed energy 
(see Chapter 5 of [lij where T{r, r') — F^g^oir, r') with 
k = 1). The difference is that we have here a new method 
to calculate the exact input kernel in terms of the fixed 
energy phase shifts through the m-function of the auxil- 
iary Schrodinger operator. 

Note that equation (fTS)) can be transformed to the 
Gclfand-Levitan form by a further transformation of T 
and JC. 



Interpolation of the m-function 

As relation ((T^) shows, the m-function must be eval- 
uated at argument values which are different from those 
offered by the input data as showed by equation ([5|). 
Therefore, we must use interpolation theories for calcu- 
lating the m-function. These interpolation theories all 
rely upon a consequence of Simon's representation of the 
m-function [l^. It states that the m-function can be 
represented essentially by a Laplace transform. For both 
the Fourier and the Laplace transforms there are known 
interpolation formulas, and the latter, a recent work of 
Rybkin and Tuan [17], will be adopted for our purposes. 

Thus, wc shall use the following interpolation formula 
for the m-function (special case of Theorem 5 of Ref. 



i{X) iVX + ^ Cn{-iVX) ^ a„m(m(-iu^) -|- i 



n=0 



with ujjn = m + ^ and 



m=0 



(17) 



{-n)m{n + l)r 



(m 



\\2 



(18) 



(i-x) 

Cn{x) ^ {2n + l)j-, r^, a„, 

{x)n is the Pochhammer symbol: 

{x)n = x{x + l)...{x + n-l), {x)o^l. (19) 

The interpolation formula applies for locally absolute 
integrable potentials. Its domain of convergence depends 
on the specific features of the potential and we suppose 
it to be valid on the whole complex A— plane throughout 
the subsequent numerical examples. 



APPLICATIONS 

As an illustration of the above theory we show first ex- 
amples where solutions of the transformed problem are 



known analytically. Such are the cases of the constant po- 
tentials and the step potentials. Then we exhibit some 
results for cases where the solution of the transformed 
problem is not at hand, and thus the interpolation for- 
mula PT|) must be applied. The selected example for 
the latter case is the shifted truncated Coulomb potential 
problem, whose reconstruction is known to be a challeng- 
ing exercise for any quantum inverse scattering method 



Constant potentials 

Let us reconstruct first the constant potential 
jqo for r < a 



which generates the transformed potential 

Q(x) - s^a^iP -qo). (21) 

The solution of the transformed Schrodinger equation ^ 
becomes 

^(x) = Ci J,^ {V~se-^) + C2J_,^^ (^/ie-^) (22) 

from which the m-function being the logarithmic deriva- 
tive of the unique solution at the origin takes the form 



m(A) 



J 



(23) 



The (s— wave) phase shift belonging to equation ([3]) can 
be expressed from (|22p a.s {n — ^/X) 



A(k) = i tanh 



4*''s""'r(iK-i-i)Ji„(vs) 
r(i-jK)j_i„(v^) 

4»^s-''=r(iK+i)j„(V^) 
r(i-jK)./_i„(\/i) 



(24) 



Now, we are in the position to perform the calcula- 
tion at three different levels. Either we can use the exact 
m-function (j23D or the exact phase function ([M)) or the 
interpolation formula (jl7p for an approximate calculation 
of the m-function using the set of input phase shifts {Si}. 
Results of the three performances are depicted in Figure 
[T] (left panel) for two constant potentials with qg = 1.2 
and a = 1 (dashed lines), and go = 0.5 and a = 0.75 (con- 
tinuous lines). The wave number in each case is chosen 
to be fc = 1. Eleven input phase shifts (/max = 10) are 
used when the interpolation formula is applied (see Fig. 
[1]). Apparently, while there is no practical difference be- 
tween the calculations performed at the levels with exact 
m-function and A-function, the weakest point appears to 
be the use of the interpolation formula. As a comparison 
we have included in Fig. [I]as dotted line the result given 
by the HA method [Tsf. 
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where 



1.0 1.5 2.0 

r(au) 




(29) 

with the Wronskian W[f,g\{x) — f{x)g'{x) — f'{x)g{x). 
The s-wave phase function can be written as (k — \/A) 



A(k) — i tanh 



4"^, "^(m+l) „ 

-1 I ^ + r(i-iK) ^ 
_ 4"s-"'r(i/^+i) „ 
r(i— zk) 



(30) 



with 



^ J-..(7i^)W[l+, 2+] ~ J,.(7iI)W^[l+, 2_] 

J«(V^)M^[l-,2-]- J_„(ViJ)W^[l-,2+]- ^ ^ 

The results obtained at A; = 1 are depicted in Figure [1] 
(right panel) at various levels of approximation indicated. 
Also in this case, the less satisfactory result is obtained by 
use of the interpolation formula (with ^max = 20) and, as 
comparison, the result of HA method is shown by dotted 
line. 



FIG. 1. Recovery of the constant H20|l with parameters go = 
1.2, a = 1 (dashed lines) and go = 0.5, a — 0.75 (continuous 
lines) and step (|25[) potentials using exact input m-function 
(upper row), exact phase shift functions (middle row) and 
fixed energy phase shifts (lower row). The HA results are 
shown as dotted line. (For further details see text.) 



Step potential 

As a next example we reconstruct the step potential 



qi for r < tq 
q{r) — { q2 for rg < r < a 
for r > a, 



(25) 



with parameters qi — 1.25, 92 — 1.125, ro = l,a = 2. 
The transformed potential looks like 



„, , \ —soe for a; < 

Qi^) = -2. , : (26) 

—Sie for X > xq, 



and the transformed Schrodinger equation ([3]) admits the 
solution 



ipix) 



(27) 

From this the m-function takes the form 



m(A) 



(28) 



Shifted truncated Coulomb potential 

We now reconstruct the shifted truncated Coulomb po- 
tential defined by 



q{r) 



— for r < a 

a — 

for r > a 



(32) 



with the parameters a — 2 and A — 1. Requiring a 
continuous logarithmic derivative of the associated wave 
functions leads to the formula for the input phase shifts 



tan 



u'i{ka) — Ciui{ka) 
v[{ka) — C'lVi (ka) 



Ci = 



kFi{kBa,ri) 



(33) 
(34) 



with kg — k'^ + A/a,ri = A/{2kB), and the usual ui,vi 
Riccati-Bessel and F; Coulomb wave functions. 

Because in this case there is no analytic result to the m- 
function or phase function A of the transformed problem 
([3]), we have the only possibility to apply the Rybkin- 
Tuan formula P?|) for calculating the m-function. The 
calculation has been performed at two different settings 
of the wave number k = 0.8 and k — 1.0. The cor- 
responding input phase shifts calculated from equation 
([55)1 are listed in Table 1 . The results exhibited in Figure 
[5] show an improvement of the procedure with increasing 
energy as more input phase shifts are involved, according 
to the semiclassical thumb rule Zmax ~ ka. However, the 
reproduction shows, in all cases, a generic departure from 
the exact result which is attributed to the shortcomings 
with use of the interpolation formula (I17|) . 
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r (au) r (au) 

(a)/max = 1, fc = 0.8 (b)/max = 3, fc = 1.0 

FIG. 2. Recovery of shifted truncated Coulomb potential H32|l 
at two different wave numbers and sets of input data indi- 
cated. 



k = 0.8 


k 


= 1.0 


I Si 


I 


Si 


-0.2991 


- 


0.3481 


1 -0.0317 


1 - 


0.0538 




2 - 


-0.0046 




3 - 


-0.0002 



TABLE I. Input phase shifts Si calculated from (|33|l for the 
shifted truncated Coulomb potential 

DISCUSSION AND CONCLUSION 

Using the Marchenko integral equation ([B]) a new fixed 
energy inverse quantum scattering method is presented 
for the recovery of potentials bounded to a finite range 
r S [0, a], beyond which it is negligible. The method 
consists of Liouville transforming the radial variable and 
the radial wave function (see eq. @), and thus cast- 
ing the problem of inverse scattering at fixed energy to 
the inverse spectral problem of a Sturm-Liouville opera- 
tor on the half line x G [oo , 0] . The possibility of doing 
this is based on the observation that the set of in- 
put phase shifts can be related to the m-function of the 
transformed problem. In [ist the Gelfand-Levitan equa- 
tion has been used to the solution of the transformed 
problem, from which one recovers the potential in the in- 
terval limr-j-o [r, a] . However one is frequently interested 
in the behavior of the potential at the origin r = 0, there- 
fore, the usage of the Marchenko equation is proposed in 
this work, because it recovers the potential (after back- 
transforming it) in the interval limr-+a[0, r]. Thus, it has 
been assumed, in view of the finite number of input data, 
that the present method recovers the potential at the ori- 
gin more accurately than that proposed in Ref . [l3| . 

By studying simple examples we have shown that it 
is, in fact, so. In the cases where we know the exact m- 
function (or phase- function) , the present procedure is ca- 
pable to recover the potential at the origin exactly while 
the method of Ref. [13| is not. 

The latter method has been formulated [13] to accom- 
modate to usage of input phase shifts by solving a mo- 
ment problem; there is no room (yet, in the present stage 



of the formalism) to offer any entry points for usage of 
exact (intermediate) input data, such as, the m-function 
and/or the A phase function. The moment problem can 
only be solved approximately even if these intermediate 
quantities were at hand. This feature and the require- 
ment of knowledge of the transformation kernel on the 
(infinite) half line leads the HA method always to pro- 
duce a worse (re)construction of the potential at the ori- 
gin. 

The present method is conducted on a way paved by 
analytical formulas derived from standard scattering and 
spectral theory. It offers explicit entry points for apply- 
ing exact intermediate input data (if known). However, 
in proportion to the moment problem of Ref. flail, the 
present technique also needs to use approximations in 
realistic situations. The determination of the m-function 
from the input data phase shifts {6i} [related to the 
wanted fixed energy potential q{r)] can be done through 
the application of the interpolation theory of the m- 
function, an up to date topic of the mathematical physics 
today. We have adopted the theory worked out by Ry- 
bkin and Tuan [17] which is an important development 
based on an earlier study [l^ . Closer inspection of their 
formula ([T7| reveals that, strictly speaking, it might not 
be applicable to our purpose: the choice ujm = m + 1/2 
in their framework is only permitted for cases when the 
norm of the potential Q is small enough and even then 
the domain of convergence of (|T7)) is such that it excludes 
the vicinity of the positive real half-line. However these 
restrictions originate from bounds on the domain of con- 
vergence of Simon's formula (overestimating the poten- 
tial by a constant) and are expected to be relaxed when 
the exponential decay of Q is taken into account. The 
fact that our numerical examples recover reasonable re- 
sults (see Figs. [T] and [5]) indicates this. At the same 
time we recognize that an important development of the 
present inverse method relies upon the more thorough 
examination of Simon's representation and the interpo- 
lation theory of the m-function. 

Another important development of the present method 
would be the inclusion of the intermediate bound states 
into the procedure. This could be done within the above 
formalism but we left this investigation to a future work. 
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